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Abstract 

We work out the relation between Chern-Simons, 2d Yang-Mills on the cylinder, and Brown- 
ian motion. We show that for the unitary, orthogonal and symplectic groups, various observables 
in Chern-Simons theory on S"^ and lens spaces are exactly given by counting the number of paths 
of a Brownian particle wandering in the fundamental Weyl chamber of the corresponding Lie 
algebra. We construct a fermionic formulation of Chern-Simons on which allows us to iden- 
tify the Brownian particles as B-model branes moving on a non-commutative two-sphere, and 
construct 1- and 2-matrix models to compute Brownian motion ensemble averages. 
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1 Introduction 



The idea that the large N Umit of quantum field theories gives closed string theories [1, 2] has 
gained much in tractability in the context of topological strings, where, as Gopakumar and 
Vafa [3] showed -in the case of the conifold geometric transition- the associated gauge theory 
is a bosonic, topological gauge theory, namely three-dimensional Chcrn-Simons theory on S"^. 
Recently, it has become more and more clear that other low-dimensional bosonic theories may 
also play a central role in string theory computations, and in particular, building on the old idea 
of Gross [4], 2d Yang-Mills [5, 6, 7]. Furthermore, these bosonic theories can often be mapped 
to statistical mechanical systems [8, 6, 9], which provides one with useful computational tools 
[10]. Also, it is by now clear that Chern-Simons and two-dimensional Yang-Mills (in its usual 
and its quantum deformed versions) are closely connected [6, 7, 11], generalizing previous work 
which related them at zero coupling [12, 13, 14] (see also [15]). 

In this paper we work out the relation between Chern-Simons, 2d Yang-Mills on the cylinder, 
and Brownian motion found in [8, 6]. We show that for unitary, orthogonal and symplectic 
groups, observables in Chern-Simons theory on are exactly given by counting the number of 
paths of a Brownian particle wandering in the fundamental Weyl chamber of the corresponding 
algebra, or equivalently by a number of non-intcrsccting movers on a line where certain boundary 
conditions are imposed. In particular, we compute the partition function, the expectation value 
of the unknot, and the expectation value of the Hopf link. Our results suggest that Brownian 
motion might give a rather general and independent way of computing knot and 3-manifold 
invariants. 

We also comment on the relation between 2d Yang-Mills on the cylinder and Brownian 
motion in an affine Weyl chamber, already worked out in [6]. This gave a connection between 2d 
Yang-Mills on the cylinder and Chern-Simons on lens spaces that made possible the computation 
of the modular transformation properties of 2d Yang-Mills on the cylinder, which for particular 
external states reduces to the quantum-deformed 2d Yang-Mills on the sphere. Thus, these can 
be translated into modular transformation properties of certain A-model amplitudes. 

We construct a fermionic formulation of Chern-Simons theory which allows us to identify 
the Brownian particles with B-model branes moving on a non-commutative sphere. 

Finally, we show how hermitian 1- and 2-matrix models can be used to compute Brownian 
motion observbles. 

This is the companion paper of [8]. Sections 2 and 3 include pedagogical introductions to 
Brownian motion and Chern-Simons theory, respectively. In sections 3, 4 and 5 we work out the 
relation between Chern-Simons and Brownian motion in various cases. In section 5 we discuss 2d 
Yang-Mills, in section 6 we work out the fermionic description of Chern-Simons and its relation 
to topological strings, and in section 7 we introduce matrix models for Brownian motion. In 
the appendices we give technical details, including a discussion of framing and a path integral 
counting of the Brownian motion paths. 

2 Random walks and Brownian motion 

Almost a hundred years ago, Einstein gave the first mathematical explanation of the phenomenon 
of 'Brownian motion', the zig-zag- like random motion characteristic of pollen grains suspended 
in water, observed by the botanist Robert Brown around 1827. Einstein's model was a discrete 
'random walk' model, that we will briefly review. 
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2.1 Random walks 



In its simplest guise, we can regard the random walk of a particle as follows. At every tick of 
the clock, the particle can jump cither right or left with equal probability. We want to compute 
the probabilty that after time t the particle will have traveled a distance x. We set t = nr, 
where n denotes the total number of steps and r the average time between one step and the 
next one; and x = ml, where m is the number of steps in the positive x-direction minus the 
number of negative steps. The probability that after n steps the particle has made m steps in 
the x-direction is then: 

(n + [n — m)\\ 2r 

where the binomial distribution comes from the fact that the particle takes a total of | (n + m) 
steps in the positive x-direction, and ^ (n — m) in the negative direction. Of course we can 
rewrite this probability in terms of x and t and take the limit of small t,1. If we do this while 
keeping 

D = f/2T (2) 
fixed, applying Stirling's formula we get for the probability density 

P(^)dar = ^S^e--'/4^* . (3) 

An important result of Einstein's computation was the computation of the diffusion constant D 
in terms of the microscopic scales I and r. 

The probability distribution (1) satisfies an interesting property, characteristic of birth pro- 
cesses [16, 17]. If we make P into an infinite matrix Pij{t), where i and j denote positions along 
the x-direction, and m the number of steps between the state i and the state j, then we can 
write the time evolution as follows: 

P'{t) = AP{t). (4) 

A is an infinite matrix whose only non-zero entries are on the diagonal, and one position above 
or below it, that is, only the matrix elements An, ^j^j+i and are non-vanishing. This 

expresses the fact that the particle moves one step at the time right or left, but cannot jump. 
The intial condition is obviously: 

P(0) = id . (5) 

Notice that, if we take r — 0, this results in the delta- function shape of the continuous case. 

Equation (4) is really the infinitesimal version of the heat equation satisfied by the continuous 
density (3). There is a nice reformulation of this equation as a difference equation [18]. Writing 
An = —(Xi + Hi), Aij = Aj if j = i + 1 and Hi \i j = i — 1, the matrix determines a system of 
polynomials by means of the recurrence relations 

Qoix) = 1 
-xQq{x) = + hq)Qq{x) + \qQi{x) 

)- {Xn+ fJ'n)Qn{x) + XnQn+l{x) ■ (6) 

Karlin and McGregor further showed that for any system of this type there exists a positive 
measure ^ for which the following orthogonal relations hold: 

/•oo r. . 

/ d'iPix)Qiix)Qj{x) = ^ (7) 

Jo TTj 
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where the moments arc given by ttq = 1, 7r„ = ^^tttF^tt^- Such a measure is called a solution of 
the Stieltjes moment problem. Extending the range of integration to all reals, one can construct 
solutions of the analogous Hamburger moment problem. 

The probability matrix is now given in terms of these polynomials by 

/•oo 

Pij{t) = 7rj dij{x)e-''^Qi{x)Qj{x) , (8) 
Jo 

and an important result of Karlin and Mc Gregor is that the determinant 

P{t-i;-{i},{j}) =det Pij{t) (9) 

is positive definite. This extends to the continuous case as well, a property that we will use later 
on. 



2.2 Vicious walkers and Brownian motion probabilities 

In this paper we consider Brownian motion of a particular type: the so-called vicious walkers 
[19]. Vicious walkers are movers whose trajectories are not allowed to intersect during the whole 
motion. Let us however consider the case of harmless movers first. These perform Brownian 
motion, but their trajectories are allowed to intersect at any time. For a single mover moving 
on a line, we compute the probability distribution of going from poing x to point y in time t, 
given a diffusion coefficient D in the medium. This is given by 

p,(x,y) = ^=i=e-MV4^*. (10) 

This Brownian motion probability has some elementary but important properties [20]. We will 
list them here for this simple case since they generalize to higher dimensions and underlie the 
matrix model description of Brownian motion and Chern-Simons. 

1) It is the unique kernel of the solutions of the heat equation with prescribed boundary condi- 
tions: 

d 

-Q^Pt{x,y) =DApt{x,y) , (11) 
which is usually denoted hy K{x,y;t). Prom now on we will set -D = ^ . 

2) It can be obtained as the continuum limit of a discrete random walk, as we reviewed in 
the previous section. 

3) It is a positive-definite quantity and symmetric under interchange of x and y. The former 
property is trivial in the one-dimensional case but it is a non-trivial result due to Karlin and 
McGregor [16, 17] for non-intersecting paths in higher dimensions. This connects the probabilis- 
tic interpretation with the Stieltjes- Wigert polynomials that naturally appear in Chern-Simons 
theory [21], as we will see. 

4) Finally, K{x,y;t), being a probability, has the extensivity property 

dniz) K{x, z- t)K{z, y; s) = K{x, y;t + s) . (12) 
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Here, the integral runs over all reals . This is sometimes called the Chapman-KolmogorofF 
equation, and is a consequence of the Markov property ('lack of memory' of the process). This 
condition also completely fixes the normalization of which in the case (10) integrates to 1. 
In the generalization we will consider later on, the normalization is still fixed by the extensivity 
property, but it is not equal to 1 since we are dealing with conditional probabilities. 

It is interesting to note that these expressions can be generalized along the lines of [20] 
to include first order derivatives in x in the heat equation. These roughly correspond to an 
external gravitational field, and the solutions are given in terms of Bessel functions. It would 
be interesting to see whether there is a gauge theory analog of such terms. 

The generalization of the above for more than one particle is obvious: 

PtA^^y) = (2j)iV/2 e"^, (13) 

and N is the number of particles. This is the case that the particles are non-interacting, and it 

can equivalently be regarded as the product of the probabilities of N single movers on a line, i.e. 
the probability for N movers on a line to start at positions yi, ■ ■ ■ ,yN and end up at xi, . . . ,xn 
after time t. 

Prom the point of view of the heat equation, however, (13) is not quite the higher-dimensional 

generalization of the heat kernel of (11), where now the Laplacian would be an A^-dimcnsional 
fiat space Laplacian. The fundamental solution is only obtained after we take the determinant, 
as we will now explain. 

Prom the point of view of the heat equation we are naturally led to consider vicious walkers, 
that is, walkers whose trajectories do not intersect at any time (see Fig. 1). The quantity 
that we are then interested in is the probability of going from a state xi,...,XAr to a state 
yi ) • • • ) Vn in time t during which the particles perform free Brownian motion but are otherwise 
non-intersecting. That is, we do not enforce the non-intersecting condition by adding a force by 
hand, but consider rather those paths among all the possible paths such that the particles do 
not intersect. Their quotient then gives the desired probability density. Further, we label the 
particles such that at time t = 0, xi > X2 > ■ ■ ■ > xn. Since the trajectories are non-intersecting 
and the particles are distinguishable, this condition remains true at all times. 

It is this condition xi > X2 > ...X]^ that allows us to think of the motion in terms of a 
single particle moving in the fundamental Weyl chamber of U{N). Indeed, if Xi now instead 
of labeling the position of a particle on the line, labels the ith coordinate of a particle in an 
A?^-dimensional space, then the motion is restricted to the fundamental Weyl chamber of U (N) . 
Later on we will develop this point of view further and see how it generalizes to other groups. 

The probability in the non-intersecting case is easy to compute thanks to the fact that the 
heat equation is linear and we can superpose its solutions. Indeed, we can use the method of 
images to write the number of non-intersecting walks as the number of free walks, from which 
we appropriately subtract "bad" walks. This is most easily done by thinking of this as motion 
of a single particle in N dimensions, in which case the "bad" walks are walks which at some 
point hit one of the walls of the fundamental Weyl chamber. This is explained in section 5. 

^In this, our conventions arc diflFcrcnt from those in [20], where for the non-compact range - only functions 
over the positive reals are considered. Further, the heat equation considered here corresponds to the case 7 = 
in [20], in which case the Bessel function is simply a cosine and the integration measure is simply d/Lt(z) = 
•y/2/7rda. Therefore, in the case considered here and with the conventions of [20] the kernel would be K{x,y;t) = 
l/\/texp{-{x^ + y^)/2t) cos{ixy/t). 
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Using the method of images (see also section 5), the probabihty is given by: 

PtM^^fJ-) = (27rtV/^ e-^^ det(e^'^^/*)i<^<j<jv 
1 

det(e 2t )i<i<j<7v , (14) 



(27rt)^/2 



where we now labeled the initial and final posititions as two vectors A and fi in M^. 

It is easy to see that this probability vanishes, as it should, when any two particles hit each 
other or, alternatively, when the particle moving in the fundamental Weyl chamber hits one of 
the walls. Indeed, in that case //,,■ = fii+i for some z, and the determinant vanishes identically. 

From the representation in the last line of (14) we immediately see that the interacting 
probability distribution inherits from the free Brownian motion case its homogeneity, that is, it 
is invariant under simultaneous constant shifts Aj — ^ Aj + c, /Xj — ^ /Xj + c. In the next section we 
will study some more properties of pt^Ar(A, /u). 

Let us briefly discuss how this interacting probability satisfies the requirements l)-4). As for 
1), consider the AT-dimensional heat equation 

^/(Ai,...,A;v) = ^A/(Ai,...,Aiv). (15) 

We are interested in antisymmetric solutions of this equation, /(Ai, . . . , Aiv), with boundary 
condition /o(A) at i = 0. The unique solution with these boundary conditions is [22] 

1 ^ r 

/(^) = ^ n / d/x^ K{X, /x; t) /o(/x) (16) 



and 



K{X, n- 1) = j^-^ det(e-(^^-'^^) Z^*) (17) 



is precisely the Brownian motion density Pt,Ni\ A*)- 
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Properties 2) and 3) were already discussed earlier, so we will skip the details. Positivity of 
the kernel K{x, y) = pt,Nix, y) was proven in [17], and its symmetry with respect to the vectors 
X, y is clear from the explicit expression. 

The extensivity property deserves a few more comments. In this case it reads 

Pt+s,N{x,y)= / zpt,N{x,z)ps,N{z,y) (18) 
Jc 

as we prove in appendix A. This simply expresses the fact that the probability of going from x to 
y in time i + s is the same as that of going from a; to y in time t, times that of going from y to z in 
time s, where we integrate over all possible intermediate states y. Notice that the normalizations 
in (14), and in particular also the fact that we integrate over the fundamental Weyl chamber C 
only and not over M^, are crucial for the proper normalization of (18). By definition, pt,N{x,y) 
-being a conditional probability- does not integrate to one, as not all possible paths from x to 
y are non-intersecting. 

Let us anticipate that (18) is equivalent to a matrix model computation of the partition 
function of Chern-Simons theory on S^. Actually it is more general than that, providing also a 
computation of Wilson loops. 

Notice that Fisher [19] includes an additional "framing" factor: 

Pt,N{>^, f^) = e-^-*^^-l^ e-^^ det(e^''^^/*)i<,<,-<Ar . (19) 

This corresponds to the total partition function, and is related to the total weights of the 
probabilities in the lattice of the underlying discrete random walk. 

In the next few sections we show how these probabilities are related to Chern-Simons quan- 
tities. 



3 Chern-Simons theory as Brownian motion: the U{N) case 
3.1 Chern-Simons theory and surgery 

In this subsection wc briefly review some aspects of Chern-Simons theory in three dimensions 
that we will need in what follows. 

Recall that Chern-Simons theory is a quantum field theory whose action is built of a Chern- 
Simons term involving as gauge field a gauge connection associated to a group G on a three- 
manifold M [23]. The action is: 

S{A) = ^ [ Tr( A AdA + l AAA AA) . (20) 
4:TT Jm ^ 3 / 

k is the coupling constant and the trace is taken in the fundamental representation. This action 
is invariant under gauge transformations; however, in order for the theory to make sense under 
large gauge transformations the parameter k needs to be an integer. In the partition function 
one then integrates over all possible gauge field configurations: 

Z = J VAe'^°3 . (21) 

At large k, the action (20) gives the equation of motion F = 0; classically, therefore, the 
Chern-Simons path integral is dominated by flat connections. 
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Notice that in order to define the action (20) there is no need to choose a metric; indeed, this 
theory is diffeomorphic invariant at the classical level, and depends only on the global properties 
of the manifold on which it is defined. This remains true at the quantum level, up to a subtlety. 
To evaluate the partition function, one needs to fix the gauge symmetry. In doing so, a choice 
of metric is required. This by itself does not present a problem; although the resulting terms 
are not topological invariant but depend on the choice of metric, one can add a counterterm 
(even though there are no divergences involved) that requires a choice of metric, such that the 
total combination is a topological invariant. In doing so, however, one needs to fix a choice 
of trivialization of the tangent bundle, in other words a framing of the manifold. Thus, the 
resulting theory is a topological invariant of framed manifolds. 

Non-trivial observables in Chern-Simons theory are Wilson loops: 

Wx{C) = TrxVexp [ A^dx^ . (22) 
Jc 

This Wilson loop is labeled by a representation A of the gauge group, and a closed loop C in 
M. This knot can have non-zero crossing number, which is the minimal number of crossings 
needed when we draw a two-dimensional picture of the knot (for an introduction to knot theory, 
see [24, 25]). Knots are topological invariants, in that they depend only on the topology of the 
manifold and on the topological class of the embedding of the loop in the manifold. Indeed, 
again one can see from (22) that the definition is independent of the choice of metric. We will 
later see how to compute Wilson loops in practice. 

It is clear that one can also consider expectation values of products of Wilson loops on knots 
that consist of several links. In that case, every link comes with its own representation Aj, and 
one should sum over the different links in (22). For knots consisting of n links, the associated 
invariant is 

Wx,...x^ = {WxACi)...WxACn)) . (23) 

This expression factorizes if the links are disjoint. The classic result of [23] is that these quantum 
field theory expectation values can actually be computed exactly, and are in fact given by 
certain polynomials (actually rational functions, when properly normalized) in the exponentiated 
coupling and rank of the gauge group, reproducing known knot invariants in the mathematical 
literature. 

Let us ilustrate this topological invariance in the abelian case. Consider for simplicity the 
case of two loops, Ci and C2. The equations of motion in the presence of these sources take the 
form: 

e^^'^d^Ax = J^' , (24) 
k 

something reminiscent of Maxwell's equations in four dimensions, where the gauge potential is 
playing the role of the magnetic field and the source for the Wilson loops is in this case 

J^{x) = ei / ds -f- 5(3) (x - y{s)) + 62 dt — 5^^^ {x - z{t)) (25) 
Jo ds Jo clt 

where ei and 62 are the charges of the particles going around the loop, and y^{s) and z^{t) are 
embeddings of the loops in M^. The above can be easily solved giving the gauge potential as an 
integral around the loop. In the classical limit we get: 
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(26) 



and filling in the value of the gauge field this gives: 



and the sum runs over a = 1, 2. x(Ci, C2) is here the Gauss linking number: 

x(Ci, C.) = ^ £ dx- £ d2/^ '-^'^^^ ■ 

As long as Ci and C2 do not intersect at any point, the Gauss linking number depends only on 
topology. It is given by: 

x(Ci,C2) = ^5^e(p), (29) 

p 

where p are the points where (on a two-dimensional picture of the link) the links cross, and 
we should not count self-intersections of the knots, which appear on projection of the knot on 
the plane. e{p) is a sign assigned to the crossing, with a for an upper pass and a —1 for 
an underpass. Notice that when the orientation of one of the loops is reversed, the overall sign 
changes. 

The expectation value (27) also contains two self-interaction terms where a = b that we have 
not considered so far. The trouble with such a term is that, even though it is well-defined and 
finite, it is not a topological invariant. One obtains a topological invariant if one provides knots 
with a framing, that is, one chooses a normal vector field along the knot, or, in other words, one 
"fattens" the knots, making them to small strips, and defines the self-linking number to be the 
linking number of boths sides of the strip. This is really a point-splitting regularization. The 
final result docs actually depend on the choice of framing; different choices of framing change the 
self-linking number by an integer, and in this way we can compare different framings. Although 
there are no divergences involved, this is analogous to a choice of counterterms in quantum field 
theory. This discussion actually extends to the non-abelian, full quantum case, and we will quote 
the result here for future reference. A Wilson loop with n links in representations Ai, . . . , A„ 
changes according to: 

(W^Ai...aJ - e2-S"=i^»'^(^>)(T^,,..,J , (30) 
where h{\) is the conformal weight of the representation A in the WZW model: 

C being the Casimir operator and g the rank of the gauge group^. Sj is an integer, denoting the 
number of units by which we change the self-linking number of the zth knot, according to the 
counting of crossings described above. 

Let us briefly recall canonical quantization of Chern-Simons [23]. This gives an alternative 
viewpoint to path integral techniques, and in particular it allows us to compute Chern-Simons 
quantities in more complicated manifolds. The idea is to construct a three-manifold M by 
surgery, that is, by cutting the manifold into tractable pieces and then gluing them together. 
Locally, the manifold looks like S x M, where we regard M as "time" , and by fixing temporal gauge 
one easily sees that the Hilbert space associated to S is the Hilbert space of flat connections on 
S modulo gauge transformations. We refer to [23] for a discussion of this. 

^For a summary of our group theory conventions, see appendix D. 
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Wc will be interested in the cases where one can obtain M by gluing together two solid tori 
after performing some diffcomorphism. The Hilbert space associated with the boundary of each 
solid torus is then the Hilbert space of the WZW model [23], which is the space of integrable 
representations of highest weight at level k. In particular, the vacuum state, corresponding to 
the solid torus without insertions of Wilson lines, is represented by the Weyl vector \p) (see 
appendix D), which is the trivial representation or an empty Young tableau. Inserting a Wilson 
loop in representation A around the non-contractible cycle of the solid torus gives a state in the 
Hilbert space, |A), where A is integrable. 

Let us consider how to get from this. To get S^, we need to identify the A-cycle of one 
of the tori with the i?-cycle of the other one, and viceversa. In order to do that, we act with 
the SL{2, Z) generators S and T: 

1 1 
1 

To' ) ■ (-) 

where T acts on the complex structure of the torus as r ^ r + 1, that is, it shifts the y4-cycle 
by a i?-cycle, and 5* acts as r —1/''"; i-e., it sends A to —B and B to A. From the latter 
operation we clearly get an after gluing both tori. Thus, we have to compute the matrix 
element 

ZcsiS^) = {p\S\p) = Soo . (33) 

That is, we need to represent S on the states |A). This is the well-known modular matrix of the 
WZW model, which we quote here together with the other generator T for future reference: 

^ wew 

2niC(X) 27vic 

Tx;, = 6xf^e^(>'+^) 24 . (34) 

The central charge is c = kdimg/{k + g), C(A) is the Casimir of the representation A, A_|_ is 
the set of positive roots, P is the weight lattice, and is the coroot lattice. 

In the gluing we could have acted with additional external factors of T. However, this only 
amounts to a change of framing of the S^, as discussed earlier. We now see that framing factors 
are proportional to T as in (34). In what follows, we construct with the transformation TST. 

By considering more general SL{2, TL) transformations, we can construct more generic man- 
ifolds such as lens spaces and Seifert manifolds. The most general manifold that we can obtain 
from gluing of two solid tori (that is, by doing surgery on an unknot inside S^) is by applying 
the SL% TL) matrix 

where — gr = 1. The representation of this matrix on afHne characters is [26, 27, 28]: 

_ (^Sgng)l^+I ^_idvr^(^(p,,)^ |P/Q|^/^ ^ ^i^^^^-^\p\^-2X{ln+wi£)^s{ln^wt,f] 

(36) 
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where Q is the root lattice, ^{U^''^^) is the Rademacher function 



p r 



p + s 



-12sip,q) , 



(37) 



and s{p,q) is the Dedekind sum 



1 

n=l 



im imp 
— cot 

q q 



(38) 



In section 5 we will work out in detail the case of the lens space S^/Zp. More general manifolds 
can be obtained by surgery on more general knots. See for example [26, 27, 28] for Seifert spaces. 

Finally, wc give the explicit expressions for the partition function and Wilson loops of Chern- 
Simons in 5"^, which wc will use later. The partition function on for gauge group U{N) is 
given by working out (33): 



-^cs('S'^ 



{k + N)^/^ 



N-l 



k + N' 



(39) 



The unnormalized expectation value of a Wilson loop around the unknot C is: 



Wx{C) 



1 



{k + iV)^/2 J 



sin 



i<3 



■n{\i - \j) 
k + N 



(40) 



3.2 C hern- Simons on S^: the probability of reunion 

In this section we work out in detail the relation between Brownian motion and Chern-Simons 
theory for the case of the partition function, which corresponds to the probability of reunion of 
random walkers. 

We will now evaluate the above probability density in the simplest case: the probablity of 
reunion, and see how this gives us the partition function of Chern-Simons on with gauge 
group U {N) . The probability of reunion is defined [19] as the probability of, given that we start 
with a configuration where all the movers are equally spaced, 

\j = {c-j)a (41) 

ending up with the same configuration after time t: jij = Xj. The special boundary condition 
(41) we will call Aoj. Notice that a is the initial spacing between two movers: 

Xoj+i - Xoj = a (42) 

and at this point it can be positive or negative. Since we arrange our coordinates such that 
Ai > A2 > . . . Ajv, wc will take a to be positive. Let us here emphasize that, even though in the 
discrete case a would naturally have an interpretation as a lattice length and hence a minimal 
distance between two movers, that is not the case for Brownian motion, where two particles can 
be arbitrarily close to each other as long as they do not intersect. Notice that to get the trivial 
representation, labeled by the Weyl vector, we have to take (see appendix D) 

N + 1 

^Oi = Pj = J (43) 
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which in particular satisfies X^^i pi = 0, & property that we will use later. However, as men- 
tioned earlier, the overall constant is irrelevant as it only corresponds to a constant shift. In 
this and the next subsection we do not keep track of coupling-independent phase factors. 
In appendix C we show that the probability of reunion is given by 



N 

N-k 



Pt,iv(Ao,Ao) = ^^^^n(l-9') 



1 ^ 

(44) 



(27r0^/2 



k=l 



where q = e If we compare this with (39) (see also (107)), we see that up to normalizations 
both expressions are the same^ if we set = 1 and 

Then we get 

Zmm = Pt,Ar(Ao, Ao) (46) 

and in canonical framing 

ZcsiS') = e-i^^^^(^'-iVt,^(Ao, Ao) . (47) 

This is the basic result in [8] concerning the relation between the partition function of Chern- 
Simons on and the probability of reunion of vicious Brownian movers. 

3.3 Wilson lines: the unknot and the Hopf link 

The next non-trivial case is to consider Brownian motion where we fix the initial positions to 
be equally spaced, but leave the final positions arbitrary. Thus we fix 

11 = pa (48) 

where p is the Weyl vector (see appendix D) and leave A arbitrary. The resulting determinant 
(14) can now be readily computed. In a first step we use the standard Vandermonde determinant 
formula (158). In a second step, we use the identity (161). The result is: 

detie-^^^-^^-y'/'') = e-^-w^(^^-i) []2sinh (^^a] . (49) 

i<j ^ ' 

Of course, we have simply proved the Weyl denominator formula for the symmetric group. As 
a simple check, filling in A = p gives us 

V^Ap. P) = e- n 2 sinh(^) , (50) 



^Notice that at this point the relation —l/t = Qs could involve a plus sign instead of a minus sign. This would 
result in a different choice of framing, which is therefore immaterial for the interpretation as the partition function 
of Chern-Simons on S^. We will see, however, that this arbitrariness is fixed by comparing our expressions with 
the WZW modular 5-matrix, which will give the minus sign. Notice, however, that ultimately this is irrelevant 
as all expressions are holomorphic in the coupling. In particular, the relation l/t = Qs should work as well if we 
identify the Brownian motion probability with the complex conjugate of S. 
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which reduces to (14) after working out the product and fihing in the value of t. 

What quantity could the more general equation (49) possibly be in Chcrn-Simons theory? 
In fact it is nothing but the expectation value of a Wilson loop on in a representation A of 
U{N), as we can see from (40). This Wilson loop winds an unknot, which in is the unique 
knot with one link component since there are no non-trivial cycles. The Wilson loop is however 
not in canonical framing. To see this, it is easiest to write 

p,,^(A,p) = eW^(^^-i)(0|T5r|A) (51) 

where the modular matrices T and S are given as in (34). In appendix E we will analyze this 
framing in detail. 

Having written down the representation (51), it is now straightforward to generalize this to 

the probability density with arbitrary initial and final states: wc simply compute the operator 
TST within arbitrary external states A and /x, corresponding to the two boundary contitions: 

p,,^r(A,/x) = eW^(^^-i)(;,|r5r|A), (52) 

as is easily checked. 



4 Brownian motion with a wall: SO{N) and Sp{N) Chern-Simons 

An interesting generalization of the Brownian movers on a line is to see what happens when we 
impose an additional restriction by inserting a wall at the origin, A = [29]. Thus we consider 
Brownian motion of N non-intersecting movers on the half- line M+, that is, the coordinates of 
the movers now satisfy Ai > A2 > . . . > AjV > 0. We immediately see that this is the Weyl 
chamber of SO and Sp, and we will return to this interpretation momentarily. 

Using the method of images, we can again compute this probability from the free case (10). 
Namely, on top of the reflections about the forbidden walls Xi = Xj, i ^ j, we need to reflect 
around the origin. We get 

N 

Pt,N{^,fi) = J2 E (-1)^^+-+^" \{p'MX,{i) - f^i) (53) 

(TESn £1...£JV=±1 i=l 

where p°(A,/x) = -^^^ q-{>^-iA'^ I'^t is the free Brownian motion probability for one particle. This 
can be rewritten as: 

PtA\i^) = jA^y^Y.^-) E (-i)-+--^-eS£i-^^«'^^/* 

^ ' aeSff £i...£jv=±l 

1 22 ^ 

= (^^2^-"^ E 6(a)n(e^^«'^'/*-e-^^«-^/*) 

1 det(e^'^^/* - e-^'^^/*)i<.<,<A. . (54) 



(27r0^/2 



2t 



We now come to the connection with the SO{N) and Sp{N) groups. (53) is really a sum 
over the Weyl chamber of -Bat. Notice that the Weyl chambers of _Bjv {S0(2N + 1)) and Cjv 
{Sp{2N)) are the same, and indeed we can interpret (54) as Brownian motion probabilities of a 
single mover moving from A to /Li in the Weyl chamber of SO(2N + 1) or Sp{2N). 
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Asymptotic expressions for (54) were given in [29] at large t for the case of initial equal 
spacing condition. However, wc can compute (54) rather explicitly using the Wcyl denominator 
formula. Remember that the equal spacing condition amounts to setting A to be the Weyl vector, 
X = p. Explicitly, 



N 



p = ^(AT + 1 - i)ei for Sp{2N) 

i=l 

^ 1 

p = ^(N+--i)ei ioT Spi2N) . (55) 

i=l 

Thus, psp = Pso + 5 (1; Ij • • • 5 !)• Obviously, although the Weyl chambers are the same, their 
root systems are different. In particular, some of the positive roots differ by a factor of 2. That 

gives the factor of 1/2 difference in the above formulas for the Weyl vectors. Now setting the 
initial condition X = p equal to the Weyl vector, both for SO and Sp groups we can rewrite (54) 
as: 



a>0 



e^^2l^l+n«-h^. (56) 



i2ntr^ 11 2t 

We now have to specify explicitly what the positive roots aij are. For SO{2N + 1) they are: 

{cj — ej}i<^j 

{ei} , (57) 

with i,j = l,...,N. Thus we get: 

N 

= (2^' " 2l l+Hsmh^^smh^^Hsmh-. (58) 

By working with the determinantal expression directly rather than using the Weyl formula, 
the two first sinh can be combined into a single Vandermonde term [30]. Having written the 
expression for generic p, we now specialize to p = p also, so that we can compare with the 
partition function of Chern-Simons. We get: 

PtMP, P) = e-l^lV*2lA|. n (smh A j (59) 

where f{k) is given in (183) of appendix F. In fact, the above is also the answer for Sp{2N). In 
that case, the positive roots are 

{cj — ej}i<j 
{ei + ej}i<j 

2{ei} . (60) 
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Thus we get the expression (58) with an additional factor of 2 in the last term. Again, specializing 
to /X = psp, we get (59) with f{k) as in formula (184) of appendix F. 

Let us now compare with the partition function of Chern-Simons on S^. It is equal to: 

ZcsiS') = Soo = ik + 9)-'/' n (61) 

a>o ^ + 5 

Up to the usual framing factor and roots of unity, this is precisely formula (56) with boundary 
condition n = p. In particular, equation (59) precisely agrees with the expressions in [31] for the 
SO{2N + 1) and Sp{2N) partition function of Chern-Simons on S^. The identification between 
the coupling and the level is given by (45), as before. The framing factor e~l^l is given by 
|p|2 = ^N{4N^ - 1) for SO{2N + 1), and \p\^ = lN{N + l)(2Ar + 1) for Sp{2N). 

It should now also be clear that the more general expressions pt^Nip,p) and pt^N{^,p) cor- 
respond to the expectation values of the unknot and the Hopf link, respectively. 

Having done SO{2N + 1) and Sp{2N), whose root systems are given by Bn and Cjv, re- 
spectively, it is now natural to look at SO{2N), which corresponds to D^. The Weyl chamber 
is now given by Ai > A2 > • • • Ajv, and in addition Aat-i > —\n- This means that A7V--1 > IAat] 
independently of the sign of Ajv- We can interpret this as Brownian motion of A'^ particles on 
a line with a wall at Ajv = such that, when the Nth particle crosses it, a mirror particle 
traveling along the mirror trajectory |Aiv| is emitted, such that Ajv-i is not allowed to intersect 
the trajectory of the mirror particle (see Fig. 2). Notice that such paths are counted twice, once 
from the particle that goes to negative Ajv, and once from the one that stays at Ajv > 0. 

The positive roots of SO{2N) are: 

af. = ei±ej (62) 

with I < i < j < N. The Weyl vector is: 

N 

p = Y,{N-i)ei. (63) 

i=l 

The probability is now easy to work out and goes as before: 

Pt,N{\ M) = ^ ^^^^e-^^ [det(e^''^^/* + e"^^/^^/*) + det(e^^'^^/* - e^*'^^/*)] . (64) 

Notice that the last term is precisely the probability for B]\j and Cm- This term vanishes if 
some Pi = 0. By construction this was never the case for and Cat, since the Weyl chamber 
precisely satisfied the condition Ai > A2 > . . . > Aat > 0. However, for this is possible for 
the iVth particle, which can indeed reach Ajv = as depicted in Fig. 2. In fact, from (63) we 
see that the Weyl vector does have a zero. Thus in computing values of pt^Ar(A, p) with either A 
or p equal to p, the second term drops out. We can now expand the determinant iox p = p to 
get: 

Pt,N{\,p) = j^—y^e 2t J__[4smh — ^smh — (65) 

which is the expected result for the expectation value of a Wilson loop for the group SO{2N). 
Filling in also X = p, what we get is the expression for the partition function of SO{2N) Chern- 
Simons on (61), as expected. The explicit expression is given in appendix F. 

A further interesting case to analyze would be that of particles moving between two walls. 
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mirror particle 




Figure 2: Two-dimensional representation of motion in the Weyl chamber of SO{2N + 1) and Sp{2N). 

5 Lie algebra wanderers: general formulation 

We can generalize the above to any Weyl group associated to a root system [32]. In particular, 
the root system can be finite (A-G series) or afHne. The affine case was worked out in detail in 
[6] and it was shown to be equivalent to 2dYM on the cylinder. 

Thus wc consider Brownian motion of a particle in the Weyl chamber C of some root system 
A. The Weyl group, W, generates M*" from the action on C, and r is the rank. Let A and 
be two points in W (which, as usual, in the finite case will be chosen such as to correspond 
to representations of the corresponding group, and in the affine case will correspond to the 
integrable representations only). We can compute the conditional probability of going from A to 
fi in time t without leaving the fundamental Weyl chamber by the method of images, as follows. 
The argument is basically a continuous version of a well-known argument by Zeilberger (see [32] 
and references therein). We denote the total number of walks^ from A to /x in time t by (A, n). 
Notice that p1{X,iJ.) is the free Brownian motion probability. 



where r is the rank of the gauge group. These are walks in W. To count the number of walks 

that stay within C, wc divide all the walks in p^{X,iJ,) into good and bad waks, according to 
whether they do or do not remain within the Weyl chamber. We will call these pt^r{^, fJ') and 
qt{X, n), respectively, and we are obviously interested in computing pt,r{X, fi). It is clear that we 
have: 



1 



g|A-M|V2t ^ 



(66) 



(27ri)'-/2 



Pt,r(A, /i) = P°,.(A, /x) - qt,ri\ /") • 



(67) 



The key step is now to use the method of images to show that the bad walks satisfy 




^This is actually a probability, but in this scctiou we will continue to talk of the number of paths; of course, 
the difference is just divinding out an infinite factor. 
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That is, for each bad walk from w{X) to // we can construct a mirror walk by reflecting with 
respect to the wall that was crossed last. Let us call the root that this wall corresponds to 
a G A. The wall is by definition perpendicular to a, (x, a) = 0. The mirror walk in question 
is then another bad walk, constructed by taking the mirror part of the walk from Waw{\) to 
the crossing point, and from there continuing to /j. along the original walk. This gives a unique 
pairing between walks in the sum (68). Since w and WaW have different sign, all such pairs 
cancel out. 

Now it is also clear that 

Pt,r{w{X),ii) = qt,r{w{X),fi) if w / id , (69) 
since w{\) always lies outside C. Combining the above three formulas, we get 

Pt{X,i^) = ^ e{w)pUw{X),f^) = -i— ^ 6(«;)e-l^-^lV2t _ ^^q) 
wew ^ ' wew 

This is the desired generalization of the non-intersecting probability, and we should note that it 
holds both if W is finite or affine [6]. 

We would next like to compute it in some natural "vacuTim" state and compare it with 
Chern-Simons theory. We will do this for the simple algebras. We take n = p where p is the 
Weyl vector labeling the trivial representation. Using the Weyl denominator formula, we get: 



/N ^ 1 |A|^+|pp T-r . (a. A) , ^, 



where a arc the positive roots. Under the identification — 1/t = 2'Ki/{k + g), this expression is 
the (unnormalized) expectation value of a Wilson loop around the unknot. 
The partition function is obtained by setting X = p: 

ZcsiS^) = -izlZ^Le-i^^-- n 2sin^ (72) 
{k + gyl^ k + g 

where we used the Freudenthal-dc Vrics strange formula |pp = ^ dimg. 

Notice that normalizing Pt,r{\p) by Pt,r{p^P) gives the quantum dimension: 

TT (73) 

We can also rewrite the above directly in terms of the WZW modular matrices. Using the 
definitions (34) and the value of the central charge c = kdmigl{k + 5), we get 

Pt,,(A,M) = eW<i-^/(T5r)A;.. (74) 

Obviously, S itself corresponds to the interaction term in the Brownian motion probability, 
and T is the Boltzmann factor, which can be seen as a propagator in the diffusive medium 
for translations over time t. It is now also clear that pt^r{_\p) itself gives the (unnormalized) 
expectation value of the Hopf link invariant with representations A and /U. 

In this language it is now much more clear why we are getting the partition function of 
Chern-Simons on S"^, and we can give an intuitive physical picture in terms of surgery that we 
will explain in more detail later. The boundary conditions A and correspond to states in the 
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Chern-Simons Hilbert space associated with two solid tori. That is, we have two sohd tori where 
the cycle that is not filled in carries a representation A, n, respectively, of the gauge group G. The 
probability tells us to propagate a random walker from A to /x but performing a transformation 
TST, that is, we glue the two solid tori with a modular transformation r — > t/(t + 1) before 
gluing. This gives us an in a framing that is not the canonical one discussed in section 3, 
but is the matrix model framing, which is the natural framing in the mirror B-model on the 
resolved conifold, as we will see later [34]. We will get back to this in section 6. We discuss 
other framings in the appendix E. 

The matrix model will be analyzed in great detail in section 7. 



5.1 Lens spaces and q-deformed 2d Yang-Mills 

We now discuss the case of the lens space /Z,p. This can be obtained by taking the afHne Lie 
algebra above, with an important subtlety that was discussed in [6]. 
The afHne Weyl chamber is 

W = W , (75) 

where T denotes translations in the coroot lattice . The Brownian motion probability in this 
fundamental Weyl chamber is then 

1tA>^^t') = E E ^He-^l^+^-""'l' . (76) 

For real positive i, this is well-known to converge; that can for example be shown by rewriting 
it in terms of 0-functions [6]. When 1/t is times a rational number, however, as is the 
case if we want to compare with Chcrn-Simons theory, the sum is not well-defined because the 
exponential is periodic. As explained in [6], a finite expression is obtained if one mods out by 
this periodicity. Alternatively, one may regularize — 1/t = 2'Kip/l + e, where p and / are integers 
and e real positive, and factor out an overall divergence, which goes like e"''/^. This is useful if 
one wants to compare Chern-Simons on lens spaces and 2d Yang-Mills on the cylinder. Indeed, 
the appropriately normalized (76) gives the partition function of two-dimensional Yang-Mills 
theory on the cylinder with the corresponding gauge group, where A and ji label the holonomies 
of the gauge field at the two ends of the cylinder. More precisely [6]: 

^2dYM(5, 9'\t)= (-yi^/Q""! g^^^(A, ^j) , (77) 

where we used the fact that the partition function is a class function of g and g' to conjugate 
them into and , respectively, and we normalized A' = 27rA/i, = 2niiJ,/l. The above was 
used to find the modular transformation properties of the partition function of 2d Yang-Mills 
on the cylinder [6]. It was further shown that: 

<f;|^^(^)'^.vMto,9';0. (78) 

up to a phase, and t' = {27ri/l)'^t. This generalized previous relations between 2dYM and 
WZW matrices [12, 13], to the case where the coupling times the area is 27ri times a rational 
number. As discussed in [6], the different normalizations are just what one expects, and the 
overall divergence is obviously related to the fact that the Hilbert space is infinite, while we are 
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imposing the periodicity mentioned before. This is however only an overall divergence that can 
be renormalized. 

One could however ask whether it makes sense to keep t real and still talk about modular 
transformation matrices. In fact, in [7] it was shown that this is possible, and defines a quantum 
deformed 2d Yang-Mills which is defined for values of gg which are not 27ri times a rational 
number^. As shown in [7], in that case the Hilbert space is infinite and one then typically sums 
over all representations and not only the integrable ones. The theory in question is the BF- 
theory introduced in [13], where the scalar cf) (the B of BF) is compact. The partition function 
was shown to be: 

A 

In [7] this was called quantum deformed 2dYM because S'ao are the quantum dimensions rather 
than the ordinary dimensions which appear in 2dYM. Let us consider the case of the sphere, 
g = 0, with 6 = 0. In this case, this is just a particular case of the partition function of ordinary 
2dYM on the cylinder where X = n = p, and so the results of [6] relate this to qt,r{PiP)i (76). 
Such an expression was also found in [7]. In [6] it was also noticed that 2dYM on the cylinder 
and the BF theory of Blau and Thompson are related, and from the above we see that indeed 

■^2dYM(P,P) = -^qYM('S'^) ■ (80) 

Notice that this holds for any value of p. The modular transformation properties of this theory 
were worked out in [6] and we can now obtain from them the modular transformation properties 
of the corresponding A-model topological string amplitudes. It would be worth to study this in 
more detail. Notice also that an explicit expression for .^2dYM(A, /x) with its modular transfor- 
mation properties are also available [6] . It would be interesting to see if these can be interpreted 
as insertions on the sphere on the q-deformed 2dYM side. 

Finally, let us remark that it is possible to obtain matrix model expressions for quantum 
deformed 2d Yang-Mills, which are useful to compute the large N expansion. In some simple 
cases, it is not hard to sec that these are the same Chern-Simons matrix model expressions, 
analytically continued in the coupling. We will get back to this issue in the near future [35]. 

6 Brownian particles as mirror B-model branes 
6.1 The open topological string side 

Chern-Simons on describes A-model open topological strings on the deformed conifold T*S^: 

xu + yv = p. . (81) 

In this section we briefly recall some aspects of the construction that we will use later on (see, 
for example, [34]). 

It will be useful to recall the canonical quantization of Chern-Simons on S'^ [36, 34]. The 
Hilbert space is as usual the Hilbert space of the torus. Denote the holonomies of the gauge field 
around the A- and S-cycles by u and f , respectively. On the torus one can always conjugate 
the gauge field into the Cartan subalgebra [36], and so 

[uhVj\=gs6ij (82) 

^This theory computes A-model amplitudes on a Calabi-Yau which is a sum of line bundles over a Riemann 
surface, where the Riemann surface is interpreted as the Riemann surface where the quantum deformed 2d Yang- 
Mills lives. 
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where i, j run from 1 to N, and as usual gs = 2m/{k + N), thus incorporating the shift 
k ^ k + N . On the soHd torus, however, as explained in [34], large gauge transformations that 
shift V in the coroot lattice are not a symmetry, and so u and v are non-compact variables. The 
partition function is then given by 

^cs(53) = (0.|e^"'/2^»|0,) , (83) 

where |0t,) is the vacuum state, i.e. a solid torus with a trivial representation. Writing the above 
in the «-basis then gives the matrix model expression for the partition function: 

/" dm . . . duN eS"=i -Q gi^i^2 (84) 

up to normalization. 

It is useful to consider the mirror of this. It is given by resolving the conifold geometry 

a;y= (e"-l)(e^-l) (85) 

on the open string B-model side by blowing np a at x = y = u = v = [37]. One then 
has N Dl-branes wrapping the P^, whose action is given by the dimensional reduction of the 
holomorphic Chern-Simons action to the world-volume of the Dl's [38, 39] which describes the 
normal directions to the branes in the Calabi-Yau^. As shown by Dijkgraaf and Vafa, this theory 
reduces to a matrix model, whose solution gives the Riemann surface appearing in the deformed 
closed string geometry [39]. This is well-known and we will not review it here. 

More generally, if we have some gluing operator U^^'^^ G SL{2, Z) which is built out the 5" and 
T modular matrices, one again obtains a matrix model on the B-model side. In Chern-Simons 
theory, the general expression for U is known and it was given in (36). 



6.2 Brownian motion fermions and the B-model 

The Brownian motion problem is an intrinsically dissipative one. It is therefore natural to ask 
what is the quantum mechanical system associated to it, that is, the quantum mechanical system 
whose amplitudes are given by the Brownian motion probabilities. Since the probabilities are 
positive quantities and the amplitudes are complex numbers, wc do expect that some analytic 
continuation will be involved. In this section we construct such a quantum mechanical system, 
and we find that it is given by a system of fermions on a non-commutative space, which can be 
identified with B-model branes [41, 42]. 

Our aim is to set up a quantum mechanical problem such that 

Pt{\,iJ,) = {fi\U{t)\\) . (86) 

We would like to find the states |A), and the operator U{t) which gives the time evolution. 
Keeping in mind the role of A and fi as holonomies of the gauge field around the two cycles of 
the torus and the canonical quantization of Chern-Simons reviewed in the previous section, it is 
natural to require 

[>^i, fJ-j] = 9sSij (87) 

but to still regard Aj and fii as the coordinates of a system of N particles, like in the original 
Brownian motion problem, since wc arc interested in describing motion of quantum mechanical 

®For a careful treatment of such non- holomorphic embeddings, see [40]. 
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particles. Furthermore, recalling the non-intersecting condition, it is now natural to take the 
particles to be fermions rather than bosons. That is, if we have a multi-particle state of the 
form 

V'i(Ai)¥^(A2)---9?(Ajv) , (88) 

we consider Slater wave-functions: 

= -^^^ det{(pi{\j))i<i<j<N . (89) 

In what follows we will also use the short-hand notation 

|A) = det(|A,),) . (90) 

Of course, we should check that this is an orthonormal set of states of antisymmetric wave- 
functions. This is easy to check: 

{X\X' ) = det 5{Xi- X'j) . (91) 

-1/ d^A|A)(A|A') = |A') (92) 
iV! Jt^n 

so really 

^ [ d^A|A)(A| = l. (93) 

The factor of A^! comes from the natural measure induced by the integration over the Weyl 
chamber: 

/ d^A/(A) = -i^ / d^A/(A) (94) 

where C is the set oo>Ai>A2>...>A7v>— oo, and / is a symmetric function. 

Now we want to be able to represent (87) in this fermionic Hilbert space. This is of course 
just a Fourier transformation. For a single particle we have 



We also have 



^iifi) = ^== dA e^'^/^'^iX) . (95) 

Obviously, we then have 

= \n/2 I d^Ae^W..^(A) = 1. ^ f d^A det(e^''/^0^(A) • (96) 
We immediately compute 

Fourier transforming back gives: 

V'(A) = (-27r/5,)^/2 f d^^e-^'^f'/ss^^i^^ . (gg) 



We are now in a position to compute the operator U{t). By definition, it will satisfy 

U{t + s) = U{t)U{s) (99) 



22 



or in other words 

(/x|C/(i + s)\X) = [ d^z. {fi\U{t)\u){u\U{s)\X) . (100) 
Jc 

Computing {iJ,\U{gs)\X), we find 

(Apfe)|,)^^-4gi 
Since this holds for any state, it also holds as an operator equation that 

Of course, we could just as well define new states 

^(A) = e-^'/2^^(^(A) . (103) 

These transform as 

(/)(A) = (-271/^3)^/2 /" d^Me-^'/2^^e-^^'^/«^0(/x) , (104) 

and their inner product gives the Brownian motion probabiHty. 

Prom (97) we see that the matrix elements (A|/x) again give the 5L(2,Z) 5-matrix. This is 
a non-trivial fact and justifies the choice of fermionic wavefunctions. On the other hand, U is 
a product of T-matrices Tx\T^n. The latter is by construction so, as we demanded pt{X,n) = 
{TST)xf,. 

It is now natural to identity the wave-functions <^(A) and <^(/x), which transform according to 
(96), as the B-model fermions of [41]. In fact, in view of the discussion in the previous subsection, 
A and /j, should really be regarded as as holomorphic coordinates rather than real coordinates. 
This leads us to consider Brownian motion in the complexification of U{N), that is, to replace 
the An-i Weyl chamber by its complexified version gl^. That is, we are considering GL{N) 
rather than U (N) . In that case, u and v really describe coordinates on a non-commutative , 
and the Brownian particles behave like B-model fermions when going from one patch to the 
other. Indeed, they transform in a metaplectic representation of SL(2,Z), and gluing one patch 
to the other with SL(2, Z) transformations is done via Fourier transformation. This is true both 
for ip{X), which transforms with S, and for the redefined fermions <f{X), which transform with 
TST. In a way, this closes the chain of dualities, since we started with Chern-Simons theory 
and wc end up with the fermions describing branes of the Kodaira-Spencer theory [41, 42]. 
Notice that to achieve this we chose a different representation of the Hilbert space from the one 
considered in the previous section, where we discussed the A-model. 

The fact that A and /j, correspond to the holonomies of the gauge field around the cycles of 
the torus suggests that we should interpret the original Brownian motion description as some 
Brownian motion in field space. It would be interesting to see how this comes about in the path 
integral formulation. 

7 Matrix models for Chern-Simons and Brownian motion 

7.1 One-matrix model and the Stieltjes-Wigert polynomials 

The matrix model for the partition function of Chern-Simons on with gauge group U (N) is 
[28]: 

2-=^][i£|?»--'^'-n(^-^^) ■ w 
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The matrix model of Chern-Simons in was solved in [21] by the method of orthogonal poly- 
nomials. The key step was to realize that a coordinate transformation 

yi = log Xi (106) 

brings it into a form for which the associated orthogonal polynomials are explicitly known. With 
the transformation (106), the Gaussian potential in (105) becomes a log-normal weight function 
w{x) = e~'°s ^'/^^ [43]. On the other hand, the sinh Vandermonde repulsion term, characteristic 
of unitary matrix models, becomes the usual Vandermonde interaction Xi — xj. The relation 
between the perturbation theory of the Chern-Simons matrix model and topological invariants 
computed in a Hermitian matrix model has been discussed in detail in [28] . 

Notice that the quantum dimensions naturally appear in the expressions for the Stieltjes- 
Wigert orthogonal polynomials [43]. Of course, this translates into the quantum dimensions 
that typically appear in expectation values in Chern-Simons theory. 

With these ingredients, the matrix model is easily solved and gives [21] 

^ k=l 

where as usual q = e^" . We should note here that the relation to the partition function in 
canonical framing is^ [28]: 

= (fc + iV)^/2 n (2sin^^J 

where in the last line we used the standard relation between the topological string coupling and 
the Chern-Simons coupling constant (45). 

In what follows we will make more comments on this matrix model -which is the same as 
the one for Brownian motion- and the computation of more general expectation values. 

Let us recall that the most characteristic property of the Chern-Simons matrix model is that 
it is a unitary matrix model with a Vandermonde interaction that is a sinh rather than the 
usual polynomial Vandermonde determinant. We will explore some interesting properties of this 
type of matrix model elsewhere [35]. In a sense, this matrix model is much simpler than the 
ordinary ones because, in fact, all the integrals involved in (105) arc Gaussian after expanding 
the factors of sinh, and in particular contain no polynomial terms, which would result in the 
usual logarithmic interactions. Thus they can actually be solved trivially without recourse to 
matrix models. The only non-trivial part is to keep track of the combinatorics, which can for 
example be done using the formulas in appendix A. Physically, the interesting part of this 
elementary exercise is that it gives us a probabilistic interpretation of the partition function of 
Chern-Simons on S^. Indeed, in the appendix we prove that the Hop! link satisfies (18). Now 
recalling the identifications (47) and (45), we see immediately that, if we set x = y = p and 
s = t, the left-hand side of (18) gives the value of the partition function of Chern-Simons on 
(with the appropriate framing as follows from (47)) as a function of the coupling, whereas the 

'^Note that, in contrast with the conventions in [8], we have included a factor of e~'^'^ in the definition of 
the Chern-Simons partition function in (108). Also, in [8] we dropped overall sign factors. 
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(108) 



right-hand side gives the matrix model expression: 



(109) 

Thus, if we set x and y to the vacuum, wc see that the cxtcnsivity property of the Brownian 
motion probability gives us a derivation of the matrix model expression of the partition function 
of Chern-Simons. It is interesting to note that this partition function comes as a function of 
gs/2 rather than gg. That is due to the fact that the left-hand side of (18) (with s = t) depends 
on 2t rather than t. This is reminiscent of a renormalization group equation. In fact, there is 
no need to set the external states to the vacuum. For generic external states x and y, (18) can 
be rewritten in terms of expectation values of Wilson lines in the corresponding representations. 
Consider the more general formula (154) derived in the appendix (with /? = 0): 

/ d^ze-tl^l'pt,iv(x,z)p,,7v(^,y) = e-^(fl^l'+^l^l')pA,iv(x,y) (110) 
Jc 

where A = ast + s + t. Rewriting this in terms of (unnormalized) Wilson loops (dropping phase 
factors) 

Pt,N{>^,fi)=e-Tt(\M'+M')Wx^(t) (111) 
where we wrote the Hopf link as a function of t rather than the level, we get 



_J^^ a^(t-l)+^ |^|2_^ at(^-l)+t |^|2) 



2A 

c 



Wx;,{A) = j^d^'ve-^M' Wx^{t)W,^{s) . (112) 



It would be interesting to find a field theory interpretation of this equation. 

Let us add that, although using matrix models and the Stieltjes-Wigert polynomials to 
compute integrals of the type (105) might seem a little overdone, the power of these techniques 
should come to full use when computing more complicated obscrvables, like those involving 
general torus links. We hope to come back to this in the near future. 

The probability of survival for N Brownian movers was defined in [19] as the probability 
that the N movers survive after an amount t of time, or, in other words, that given some intial 
condition /i, after time t the trajectories of none of the movers crossed. Given that pt^N{X,^i) 
is the probability of moving from to A without intersecting, to get the probability that the 
movers have not intersected at t we have to integrate over A: 

Psurvival(A) = / d^Apt,Ar(A, /x) 
Jc 

= ^ [ d^AK^(A,^)| 
iV! J^N 

where as usual the integral is over the Weyl chamber. In the second line we made use of the fact 
that pt^N is antisymmetric in A. This can be approximated for large t as follows [29]. We choose 
t large and rescale A such that |A|/t is small but keeping |Ap/t large, and take the boundary 
condition so that l/xp/i is small (and the particles are equally spaced). Keeping only linear 
terms, we get 
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It is clear that the factor of came from the Vandermonde interaction. It is amusing 

to note that this combines with the first term to give a factor of t^^ The above integral can 
be readily done using Selberg's integral [44] and is given in terms of F-functions. After doing the 
integral, the total t-dependence is ^-^(^-i)/^ p9]. Expression (114) was also found by Fisher 
[19]. 

Notice that from the above, for large t, we have 

Pt,iv(A, m) ~ e-^^ A(A) . (115) 
This expansion can be used to compute more general observables: 

{f) = ^ f d^A|p,,^(A,M)|/(A) (116) 
where we assumed that / is symmetric in fj,. One can again approximate this for large t as 

- mi2^'~''"'-"" n ^ I d^Ae-^ n 1^. - ^.1 m . (117) 

Again, for a large class of symmetric functions one can make use of Selberg's integral to evaluate 
this explicitly. 

For antisymmetric functions, we write 

/(A) = A(A)/(A) (118) 

where / is symmetric. Then we simply get: 

</W> - M(dp^'-"<"'-'"^n^ /^^d-'Ae-^AfAlVCA) (119) 

that is the Hermitian matrix model with (3 = 2. This is the usual case, and one can use the 
results of for example [45] to compute the kernel. Its short-distance behavior is well-known: 

K{\^^) = '^^^^^^^^^^ (120) 
7r(A - //) 

for \ — ^ <^\. At this point one can make use of the usual matrix model techniques to compute 
Brownian motion averages in the approximation of large i, which corresponds to weak coupling 
on the field theory side. 



7.2 Brownian motion averages and the Hermitian two-matrix model 

We have seen how Brownian motion is related to one-matrix models. The probability pt,jv(A,/x) 
can however be interpreted more generally as an integration kernel in the context of two-matrix 
models. For simplicity we will do this for the case of U{N). Brownian motion two-matrix models 
have also appeared in [46]. 

We start with a simple example for the case N = 2. Assume we are interested in the 
average of the final distance between the two particles, (/ii — ^2) after time t, given some initial 
configuration Ai — A2. As usual, we take Ai > A2, //i > /X2- We can compute this as follows: 

ifJ-i - IJ'2) = I d/iid/X2Pt(/x, A)(/xi - 112) (121) 

^ JlJ,l>IJ.2 
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and in this case, explicitly, 



P'^^'^^ = 2^t - e^^T^) . (122) 

The normalization is given by the probability of survival 

p(/x,A). (123) 



/ 



'Ml>/"2 

So we get 

(Ai - A2) = ^ / p(m, A) (mi - //2) = Ai - A2 , (124) 

and Z was computed asymptotically in (114). In what follows we will drop such normalizations. 

One can more generally consider averages of functions in both variables, /(A, n), correspond- 
ing to a process from A to in time t, and one may also be interested in integrating over both 
boundary conditions, initial and final ones. Consider now for example: 

; (m -/^2)(Al - A2) , 3 , . 

^ (TTWW^ = ^ ^ 

where we set t = 1. From the above it is already clear that such integrals become more and 
more complicated as N becomes large, and also one needs to make sure that they converge. In 
the above, the presence of the denominators ensured convergence. Notice that in this case —but 
it holds more generally- after performing the /^-integrals all the Gaussian factors disappeared, 
and so in general one has to consider rational functions rather than polynomials. 

In the following we will show how two-matrix models can be useful in this type of computa- 
tions, by mapping the density pt(A, /i) to an integration kernel of a two-matrix model. 

Thus we consider quantities of the following type: 

N 

/ n dAfed/Xfe det{e-^^*-^^^"/^')i<i^j<N /(A, m) , (126) 

and the integration region C is as usual the fundamental Weyl chamber. The hat on / is to 
indicate that it is an antisymmetric function under interchange of any two Aj or /ij. Notice 
that this is not necessary, as antisymmetric functions do not vanish when integrated over the 
fundamental Weyl chamber (in contradistinction to their vanishing when we integrate over the 
whole of M.^). Nevertheless they are useful in that they can be extended to integrals over M.^ 
due to the identity: 

y^ndA./(A)=^y^^ndA./(A), (m) 

where C is the fundamental Weyl chamber, for a symmetric total integrand /, whereas for 
antisymmetric functions: 

N N 

n dA. /(A) = ^ E / n dA. /(A) (128) 

"'C'i=i wew "'"'(C) i=i 

where w{C) is the Weyl chamber reached by acting with an element of the Weyl group (in 
this case a permutation) on the fundamental Weyl chamber. Of course, we can always turn 
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an antisymmetric function into a symmetric one by taking its absolute value, and in that case 
one can extend the integral to R^, as we saw in the case of the one-matrix model. For the 
time being we will limit ourselves to integrals over symmetric functions, that is, antisymmetric 
functions / in (126). For these ovserbables, the ordering of the particles is crucial. Without loss 
of generality we can write 

f{X,^i)=A{X)A{^^)f{X,^^) (129) 
where A (A) is the usual Vandermonde determinant 

A(A) = n(Ai-A,). (130) 

i<j 

Thus we have 

N 

if) = N\ ndAid/x,A(A)A(M)det(e-(^-^^)'/2*)/(A,M) 



1 r ^ 

n'\ Ln n 



dXid^Xi A{X)A{fi) det(e-(^*-'^^)V2*) /(A, ^) . (131) 



Remember that /(A, /Lt) is a symmetric function. In the last expression, any antisymmetric piece 
would not contribute. We will further assume that / is rational in A, /x and chosen such that 
the integral converges. We included the factor of AT! and dropped a factor of (27rt)^/2 (''^'^ich we 
will nevertheless take into account later on) to simplify normalizations. 

Before moving to the two-matrix model, it is convenient to rewrite the above in an alternative 
way. We note the following formula for arbitrary symmetric function: 

/ d^yA(y)det(e-^^Ol<^<.■<iv/(y) = iV! / d^y A(y) eSf=i-^2/^ /(y) . (132) 

This is easily proven by writing out the definition of the determinant as a summation over the 
symmetric group, interchanging the order of summation and integration, and relabeling the y- 
coordinates. Similar manipulations are illustrated in appendix A. This results in a summation 
over Sn that gives the factor of A^!. Of course, the right-hand side is still antisymmetric in the 
property that will be crucial in what follows. 
We can now rewrite expectation values of antisymmetric functions / as 

(/) = / r[dAid/x,A(A)A(/x)e-^£i(^'-'^')V2tj(A,^) . (133) 

Again, any antisymmetric piece in / drops out after symmetrization. 

It is not hard to recognize in (131) a two-matrix model. The form (131) is the one used in 
[22] which results from the Itzykson-Zuber integral, the simpler one (133) is the one advocated 
in [47]. It is a two-matrix model interacting via a potential "^^Li XiUi, as we will now see. 

Consider then the Hermitian two-matrix ensemble [22]: 

Z = l [dL] [dM] e-nL)-V{M)+pTriLM) ^-^34) 

Averages in this ensemble are integrals of the form 

{f{L, M)) = j [dL\ [dM\ e-nL)-v{M)+pTr{LM)j^^^ ^) _ ^^35) 
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Typically, f{L,M) is a rational function of the moments TrL™, TrL". 

Wc will clearly be interested in quadratic potentials^ V{M) = M^/2i and P = 1/t, in which 
case the above takes the form: 

(/) = J [dL][dM] e-^(^-^)'/2*/(i^, M) , (136) 

for some function f{L, M). For now, however, we will leave the parameters generic. The measure 
[dL] has been worked out in [22, 47] and it contains two parts: an integral over the diagonal 
entries of L, and an integral over the angular part. Indeed, we can diagonalize L by a unitary 
matrix U as 

L = UAUK (137) 

where A = diag(Ai, . . . , A^v). This, however, cannot be done simultaneously for L and M. Thus 
there remains an integral over unitary matrices that we have to perform explicitly. We have [47] 

N 

dL = dAdU A2 (A) , dA = dAj . (138) 

i=l 

References [22, 47] then prove that the above reduces to 

= ^7?S^ / n d^^d/x. A(A) A(;.)e-^W-^('^)/(A, /.)/(A, ^^) , (139) 

(111 p^-r J i=i 

and we are assuming that /(L, M) is a function of the moments of L, M and no mixing terms. 
J(A, jj) is the Itzykson-Zuber integral, given by 

and it is the result of the integration over the angular variables. We see that up to normalization 
this is precisely the integral^ (131): 



N-l 



7 ^ ^ V A(Ai)A(/ii) ' 

Putting everything together, this gives the Brownian motion ensemble (131). 

One can use uniqueness of the solution of the heat equation with respect to the unitary 
invariant Laplacian operator on Hermitian matrices and on the Cartan subalgebra, respectively, 
to compute the coefficient appearing in the transformation from dM to an integral over the 
diagonal. What is of interest for us is the fact that the solution kernel of the former equation is 

f{L, M; t) = ^^-i^ ,-'T.(L-Mf,2t ^ (142) 



^Notice, however, that in this special case Z is not well defined: it diverges like ~ where A is a cutoff on 
the L, M-intcgrals, or, alternately, as ~ \/{fi— l/t)'^^'^ [48] as /3 — > This is however not the Brownian motion 
normalization. 

®In this formula M denotes the diagonal entries of M and not M itself. 
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and this explains the t-dependence in (140). Let us now turn to the heat equation on the Cartan 
subalgebra. This is of course the heat equation we considered before: 

^/ = iA/, (143) 

and we look for antisymmetric solutions of this equation, /(Ai, . . . , Ajv) with boundary condition 
/o(A) at t = 0. These are given by [22] 

N 

/(A) = n Uii^K{X,p.-t)Uij) (144) 

and 

^(^'/^;*) = J^2W\ det(e-(^-''^)^/^*) (145) 

is precisely the Brownian motion density pt^Ni^lj)- Of course, any symmetric solution of the 
heat equation can be obtained by rescaling the corresponding symmetric solution with the 
Vandermondc determinant. 

This is of course analogous to the case of Brownian motion in the affinc Wcyl chamber of 
?7(A'^) /two-dimensional Yang-Mills [6], where the partition function is the heat kernel of solutions 
of the heat equation for symmetric functions. 

Having found a correspondence between Brownian motion averages and one- and two-matrix 
models, we can now use the techniques available for these models, in particular large-iV tech- 
niques. In particular, we have the following correlation function: 

PM{\l,...,\N^^^) = ^|a(A) det(e-(^-«)V2*) (146) 

in terms of which we can define our Brownian motion averages for antisymmetric functions, as 
before. The normalization factor is 

Z = yd^AA(A)det(e-(^'-'^^)'/2t) _ (147) 

This is the analog of the density of eigenvalues in the presence of a coupling to an external 
source /x. It is related to the two-point kernel as follows: 

Pn{\i, . . . , Ajv) = deiK{Xi, A,) . (148) 

We supressed explicit /x-dependencc because here we are considering it as a source, but these 
expressions depend on ^. An explicit expression for K_ has been computed in [45], and the 
large- asymptotics and short distance behavior have been analyzed in [45, 49]. It is well- 
known that K can be expressed in terms of orthogonal polynomials. The generalization of this 
to the present case is in [45]. One can also integrate out some of the coordinates Aj to get a 
correlation function depending on a smaller number n, n < N, of eigenvalues. In that case the 
determinantal identity (148) continues to hold [45]. 

Notice that this kernel is not the kernel computed before; symmetric solutions of the heat 
equation can be obtained from it, and it satisfies (12). In fact, it is related to the Brownian 
motion kernel as follows: 



1 

/d^AA(A)i^(A,M) 



deiK{Xi,Xj;n) = A(A)K(A, /x^ . (149) 
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Two-matrix models are solved in [22, 47]. [50] has found exact formulas for mixed correlation 
functions of the type Tr L" Af using biorthogonal polynomials. 

Notice that the same two-matrix model was conjectured to play a role in the computation 
of correlators of branes in Calabi-Yau crystals [42] . Since in the previous section we already saw 
that Brownian motion particles can be interpreted as B-model fermions, it would be interesting 
to see if the above formulas can be used in the context of [42], and whether they have a Chern- 
Simons interpretation. 

8 Discussion and outlook 

In this paper we have shown that for a number of observables the mathematical structures of 
Chern-Simons theory and Brownian motion for non-intersecting movers are almost identical, for 
the U{N), SO{N) and Sp{N) groups and for the manifolds and S^/Zp. The reformulation of 
the classical Brownian problem as counting of quantum mechanical paths in the Weyl chamber 
of the gauge group points to a deeper connection between Chern-Simons and Brownian motion, 
namely at the level of the path integral where we count paths on field space. Indeed, from the 
canonical quantization approach it became clear that the coordinates of the particles describe 
the holonomies of the gauge field around the cycles of the torus along which one does the gluing. 

An intriguing connection with the B-topological string was found at the level of fermions. 
The natural fermionic formulation that we find for the Brownian particles is identical to the one 
in [41], thus identifying our particles with the branes of the B-model. In the case of that 
we analyzed, we found that these fermions move on a non-commutative space with two patches. 
For this we needed to complexify the Weyl chamber, thus effectively replacing U{N) by GL{N). 
Brownian propagation then translates into a coordinate transformation from one patch to the 
other which is a Fourier transformation. 

The case of lens spaces corresponds to motion in the afHne Lie algebra, and is related to 2d 
Yang-Mills theory. [6] related Chern-Simons on lens spaces to 2d Yang-Mills on the cylinder at 
a value of the coupling that is 27ri times a rational number. However, the results in [7] suggest 
to consider also the analytically continued Chern-Simons. Since integer values of the coupling 
are required in Chern-Simons in order to preserve invariance under large gauge transformations, 
the deformed theory is defined by the quantum deformed 2d Yang-Mills, which for the case of 
the sphere is the same as usual 2d Yang-Mills on the cylinder, as we showed. Thus, the modular 
transformation properties worked out in [6] should correspond to those of the A-model on the 
corresponding bundle over the two-sphere. It would be interesting to investigate this further. In 
fact, Brownian motion would seem to be more naturally associated with the quantum-deformed 
2d Yang-Mills than with ordinary Yang-Mills, due to the exponential nature of the probabilities. 
This quantum deformed 2d YM [7] is an interesting theory in its own right. The partition 
function was evaluated in [6, 7] for the sphere and the torus, and by using the manipulations 
in appendix A it should be readily extendable to the higher genus case. In particular, one can 
obtain matrix model expressions, which are useful to work out the large limit. We will come 
back to this in the future [51]. 

It is likely that the generalization of the general modular matrix U^^'^^ is simply given by the 
same expression (36), where one analitically contimics in the coupling and sums over the whole 
coroot lattice. It would be interesting to check this. In that case one can use the relation to the 
0-functions to study its modular transformation properties, as in [6], and one could study the 
reformulation of Chern-Simons theory in terms of 2d qYM for more general manifolds. 

One important issue is to extend the Brownian motion description of Chern-Simons to more 
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general manifolds and knot invariants. Particularly interesting is the case of knots consisting 
of three links or more, since so far we encountered only up to two-hnked knots. This is of 
course due to the fact that we only have two boundary conditions to our Brownian motion 
problem. This suggests that for extra links - which will carry additional representations - one 
will have to impose more boundary conditions on the problem. One could for example consider 
the conditional probability of starting with a configuration A, ending up with a configuration 
/J,, given that the intermediate situation at some time t was an intermediate state u. It would 
be interesting to check whether this associates an invariant with representations A, i/ to a 
three-link knot. It would be very interesting to work along these lines further, as this might give 
an independent and general way of computing knot invariants. It might be useful here to really 
view the Brownian motion picture as a braid diagram, and consider the associated algebra. 

An interesting open problem is to develop matrix model technology further so one is able to 
compute for example general torus links, and in particular their large N expansions. The same 
is true for three-manifolds M where M is a Seifert homology sphere, so one can work out the 
corresponding geometric transitions on T*M. As far as Brownian motion is concerned, we saw 
that in order to replace by S^/Zp one had to replace the finite fundamental Weyl chamber by 
the afHne Weyl chamber, which amounts to modding out by translations in the coroot lattice. 
It would be interesting to generalize this procedure to other manifolds, and in particular see if 
there is a general way of modifying the manifold by modifying the shape of the region where 
the particles are allowed to walk. 

Recently, the melting crystal picture was applied to Chern-Simons on directly [52]. It 
should be possible to see directly what melting crystals and Brownian motion have to do with 
each other. 

Finally, it would be interesting to understand the role of the heat equation within Chern- 
Simons theory (perhaps in terms of stochastic quantization [53], as a Fokker-Planck equation) 
and in the topological string itself. 
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A Extensivity properties of the kernel 

In this appendix we give the proof of the extensivity formula (18): 




(150) 



and work out several other identities of the same type. We first write pt,N{x,y) as: 




(151) 



so that the right-hand side of (150) is 



iV^g-E£l[(2a(i)-^i)V23+(^.'«-2/i)V2t] _ 



(152) 
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We now relabel the Zi^s and perform a transformation a a' a. The a' sum then results in a 
factor of A^!. Bringing the remaining terms together and performing the Gaussian integrals, we 
are left with: 

which is what we wanted to prove. 

In fact we can easily extend this computation to the following: 



(27r(s + 1))^/2 



c 



d^z e 2 l^l'+/?Efci ^^ pt j^{x, z) Ps,n{z, y) 

= e-2kNfl-l'+!l2'l')-2/^i:£i(-'+*-^)-^/3'*^PA,iv(x,2/) (154) 

where A = ast + s + t. We can also replace /? Xi by YliLi Pi^i the integrand, and then 
the final expression contains an additional symmetrization with respect to 

Notice that we can construct gft,r(A,/i) - or, equivalently, the partition function of Chern- 
Simons on lens spaces - from two partition functions. Indeed, up to normalizations, we 
have 

gt',r(A,/x)= (^)'^' J2 [ J'ze-'l'?+^^'^/%r{\z)ptA^,fi) (155) 

where the coupling is t' = a^t. The left-hand side is the operator ST'^S. One can also show 
extensivity of qtA\l^) itself. Analogous expressions can be derived for q-deformed 2dYM. 



B Path integral formulation of the path counting 

Computing non-intersecting Brownian motion probabilities involves picking an infinite set of 
non-intersecting paths out of the infinite set of all free Brownian motion paths. It is therefore 
natural to look for a path integral formulation, where particles are traveling from an initial state 
to a final state along Brownian paths. It is well-known that free Brownian paths correspond to 
quantum mechanical paths, in other words, the free Brownian motion probability is the quantum 
mechanical path integral of a free particle. Thus, all we have to do is take into account the fact 
that the particles are non-intersecting. By elementary manipulations, we get: 

Pt,iv(A,M) = det^"*' Vxi{t) . . .Pxiv(t)e^^[^Wl (156) 

where the determinant is taken with respect to the boundary conditions Aj and /Xj, and <S'[a;(f)] 
is the action for a free particle. It is interesting to discretize this as follows: 

Pt,N{\lj) = 

= Jim (^) 1^ dxl...dx%... 1^ dxf . . . dx^-' n' det(e^(<-4^^)^)„,(;i57) 

where we simultaneously take the limit e ^ and M ^ 0. The coordinates xf label the 
discretized coordinates, where i = I,. . . ,N and o = 1, . . . ,M. Notice that we were able to 
restrict all integration regions to the Weyl chamber, as expected. Thus, this can be seen as a 
path integral for motion inside the Weyl chamber. 

It is likely that the above expression presents an additional hidden supersymmetry when 
expressing the determinant as a fermionic integral, in the limit M — > oo, as in [54]. It would be 
interesting to work this out further. 
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C Some determinantal identities 



In this appendix we collect the steps leading to (44). We use the following identities: 

^<i-<j<N 

Yl q^^ = gi:£2(i-i)Ai . (158) 
l<i<j<N 

The first is the standard Vandermonde identity, and the latter is also easily Checked. 

We start from the last equality in (14) with Aj = /Ltj = (c — i)a and extract the Gaussian 
factors: 

det(e-^)i<,<,<^ = _J__e-E£.(^-i)^«VMet(eO-^)(-^^^ , (159) 

where we used the freedom to shift the initial and final positions to recast the determinant in a 
form in which we can apply the standard Vandermonde determinant formula (158). Then using 
both of (158), we see that the prefactors cancel and we finally get 

N 

det(e-(^-^)'«'/2*) = n (1 - (160) 

k=l 

with q = e""^/*. 

For the computation of the expectation value of a Wilson we need in addition the identity 
J](e^^ -e^O = n^^"'n2sinh^^^ . (161) 

i<j k i<j 

This is easyly checked for low values of N, and proved by induction. 



D Representations and partitions 

We first introduce some Lie algebra conventions. In this appendix we deal with the case of U (N) 
and SU{N). We start with U{N). We use standard notations (see for example [55]), where Ui, 
i = 1,. . . ,N, denotes the set of fundamental weights, and are unit vectors in R^. They are 
related by: 

i . N 

j=i i=i 

The simple roots arc given by a.j = e.j — ej+i for i = 1, . . . , N — 1. From here it follows that 
|ap = 2 and so a = a^. Now we can expand a highest weight of an irrep 

TV TV 

\ = ^\iiOi = ^{£i-K)ei (163) 
i=l i=l 

where 

1 ^ 
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is the number of boxes in the Young tablaux divided by N, and 

ii = X^ + X^+l + ...Xn . (165) 

Using this definition, it is not hard to check that (163) holds. 

A Weyl chamber is a connected set in left out after we delete all hyperplanes orthogonal 
to the roots. It is defined by 

Cy, = {X\{w\, ai)>0,i = l,...,r}, weW . (166) 

The fundamental Weyl chamber corresponds to the identity element and is denoted Cq. Obvi- 
ously, Aj > for all i = 1, . . . , if A is in the fundamental Weyl chamber. This implies that the 
£'s satisfy 4 > Defining coordinates 

N + 1 

hi = £i + - i , (167) 

we see that the fundamental Weyl chamber is the domain hi > h2 > ...hN- These are the 
variables that are usually used to label representations, and which we have labeled by A and jU 
(and we will continue to do so). They can be written as 

h = i + p, (168) 

where p is the Weyl vector 

^ = E-^ = E(^-0^^ (169) 



i=l 



and has norm 



H2 = liV(iV2-l). (170) 



The fundamental weights satisfy: 

iuuuj) = F,, = = {C-')ij (171) 



lap 



where in the last line we specialized to U (N) . C is the Cartan matrix. The inner product for 
weights is 

(A,/x) = ^Ai/x,Fy . (172) 

For U{N), the Casimir becomes especially simple in the new basis: 
C2{X) = (A,A + 2p) 

^ 1 

= \£ + p\^-\p\^ = \h\'-\p\' = Y, - - N{N' - 1) . (173) 



The dimension of a representation is given by 



dimi?= TT ^^^-^ (174) 

l<i<j<N •' 
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E A note on framing 



In previous sections we related Brownian motion and Chern-Simons quantities. Framing in 
Chern-Simons theory was already explained in section 3. Both in the case of the partition 
function and for the unknot we saw that the expressions resulting from Brownian motion were 
not quite in canonical framing, but in matrix model framing. That is, every state came multiplied 
with a Boltzmann factor. We will now explain how such factors can be modified by changing the 
boundary conditions of the Brownian motion. We will consider translations and dilatations of 
the boundary conditions. Besides the possible interpretation as framing factors in Chern-Simons, 
these properties of the probability densities are interesting in their own right. 

We start with an example that will however cover all of the physical effects we want to 
analyze. We take initial and final states that are equally spaced, but are not identical: 

Xi = (c — i)a 

fj,i = (c' — i)a' . (175) 

Let us further introduce the notation 



a = ca — c'a' 



(3 = a-a! . (176) 

We are thus set to compute the determinant (14) for the above values of A,//. This is not hard 
to do and the procedure is similar to the way we found the probability of reunion. Basically, we 
expand the exponential in such a way that we can use formula (160). The result is 

N-l 

det(e-^^^-''^)'/2*) = (>[-kP^N(N+mN+l)-Na(a+l3(N+l))\ ^ -|-|- _ ^-kaa'/t^N-k _ (^77^ 

k=l 

We see that a and f3 can be chosen such that the exponential equals a framing factor 
'^i2{k+N)^ n for some integer n. So in principle we can reach any framing we want by translating 
the initial or the final positions of all the particles by a constant factor, or by rescaling them. Let 
us however stress that these are not symmetries: the probability densities change by exponential 
factors. However, once we have computed the probability for a specific configuration, wc can 
compute that of any other configuration connected to it by translations or dilatations with the 
above transformation rules. 

Not all of the above transformations have a natural interpretation in terms of framing though; 
the translations do not. To see this, it is best to look not at the "probability of reunion" but 
more generally to the transformation of pf.Af (A, fi) under dilatations and translations. Let us first 
consider the latter. We already know that the probability density is invariant under translations 
by equal factors, A — A + c, — +c. We now consider the following: 

p*,^(A + ca,M + c'a) =e-^"=i(^'-'^')(^-<=')f-^(-'=')'l^Pt,jv(A,M) ■ (178) 

Clearly, the first factor does not have much chance of being interpreted as a framing factor 

except in an SU{N) theory where it vanishes. For U{N), the framing factors go with the second 
Casimir of the representation, not with the first (see section 3). Of course, this factor vanishes 



^°We will continue to use this name, even though with the boundary conditions (176) the particles do not go 
back to their initial positions, now even up to an overall translation. 
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when |A| = and in particular for the Weyl vector, which satisfies |p| = Ylj^iPi = (see 
appendix D). But there is no reason to restrict oneself to this type of boundary conditions, so 
this does not generally reproduce the effect of framing. 

The effect of dilatations on the external states however is much more similar to framing. This 
corresponds in (176) to a = (3. Let us now again consider the most general case of arbitrary 
initial and final states. To this end it is easiest to work directly with (14). If we rescale 



A 



Xa 



we get 



Pt,N{Xa,iJ,a') = 



|Ap 2 M^„n , , 
^« det(e 



Xifijaa' /t 



)• 



(179) 



(180) 



(27ri)^/2 

Of course, this reduces to (177) for a = /3 and the probability of reunion. However, we do 
not want to change the last interaction term, which corresponds to the modular matrix S and 
does not change under framing. At this point we can either generalize the relation between t 
and the Chern-Simons coupling to include the extra factors of a, a', or we can consider rescalings 
of time as well: 

e"— « -— « det(e^'''^/*) . 



^*'^^^''^'^=(W^ 
This now corresponds, up to an overall constant, to a matrix element 

(/x|T'^5r™|A) 



(181) 



(182) 



which is what wc wanted. For this to be a framing, c? = n and a'^ = m are integers. With a 
and a' real, notice that this can only augment the framing and not diminish it. 



F Partition function of SO and S]) Chern-Simons 

F.l ,SO(27V+ 1) and 5p(2Ar) 



Let us now give j{k). For SO{2N + 1): 
f{k) 

For Sp{2N): 

f{k) 



F.2 S0{2N) 
For SO{2N): 



'1 

N-k/2 
^N-k/2-\ 



k even 
k odd . 



{j = l,...,N) 



( N 

N 



k/2-\ 
k/2 



N -k/2 + \ 
N-k/2 + 1 



kodd<N 
k even < N 
k odd > N 
k even > N . 



(183) 



(184) 
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Again, using Weyl's denominator formula this can be brought to the form 



Soo = {k+ 5)"'^/' n ^'^^ hr ^'""^ 



2N-j-i 
2t 



i<j 



N-1 2N-1 



{k + 2N -2)-^/'^ W{smhk/2ty^^^ JJ (sinh A;/2t)-^W 



1-1 



fe=l k=N 



2N-1 2N-1 



(A; + 2iV - 2)-^/2 JJ (sinhA;/2t)^(*^)/ JJ sinhA;/2i 



(185) 



fe=l k=N 



and f{k) = N — l/2k if A; is even, and N — l/2k + 1/2 if it is odd, in other words. 



fik) = [N- 



k + 1 
2 



(186) 
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